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Abstract 

We prove that the dynamical system charaterized by the Hamiltonian 
H = XN J2f Pj + X/^fc (PjPk) 1 ' {cos[v(qj — qu)}} proposed and studied by Calogero 
[1,2] is equivalent to a system of non-interacting harmonic oscillators. We find the 
explicit form of the conserved currents which are in involution. We also find the 
action-angle variables and solve the initial value problem in simple form. 
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1 Introduction 

Motivated by the discovery [3] that the dynamical system characterized by the Hamilto- 
nian 

N 

H = J2 (PjPk)exp{-ri(qj - q k )} (1) 

3,k 

is completely integrable, Calogero [1] proposed a new Hamiltonian system with the stan- 
dard Poisson bracket and the Hamiltonian 

N-l 

H = XNP + (P 3 PkY 2 {cos[v{ q] - q k )]} (2) 

j,k=0 

N-l 

P=Y,P> (3) 

3=0 

The equations of motions of the dynamical variables are 

d i N ~ 1 

—qj = XN + /ipj — Pk* cos[is(qj - q k )\ (4) 
d i JV_1 i 

—Pj = 2fiup^ Pk* sin[z/(g j - q k )\ (5) 
at k =o 

Then by way of a lenghty anlalysis, he succeeded to show that this Hamiltonian system 
is solvable, in the sense that the above evolution equations can be solved in closed form. 
The final form of his solution take the following form: 



qj (t) = qj(0) + XNt+ 

-i arptan ( sin (ujfiNt + a - %(0)]) - sin {v[a - qj (Q)]) N m i\ , , 
V arCtan [cos WNi + a- q~m) - cos („[a - <&(0)]) + A [PA ° )l ) (6) 



Pj (t) = Pj (0) + 2-b,(0)]5 ^ cos (u\pNt + a- q,(0)]) - cos (u[a - <&(0)])J + 

2(4) 2 [l-cos(^M)] 



(7) 



Where the constants A and a are determined by the initial data as follows: 



C(0) = A cos a 




(8) 



5(0) = A sin a 



AT-l 

EPfc(0) 3 sinK?*(0))] 



(9) 



fc=0 



This part of the analysis was done mainly by solving the equation of motion for the 
functions C(t) and S(t) plus lenghty and at times cumbersome use of trigonometry. 
He then went on to study the " Behavior Near Equilibria " ( only in the case A = ) 
and by analysing this system along the lines of the standard theory of small oscillations 
, concluded that [1]: 

" .... In the case A = considered here, the general motion is characterized by a com- 
pletely periodical behavior..." and In view of the remarkable simplicity of the motion 
characterizing the model, it is natural to conjecture that there also exist a quantal version 
of this model and perhaps of some of its generalizations which is solvable and features a 
spectrum whose discrete part is equispaced....". The quantum version of this model was 
then discussed in ref . [2] . 
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This part of the analysis was done by heavy use of properties of some special matrices 
introduced for this purpose in [1] . 

The discovery of new integrable systems have always created excitement and great 
activity in mathematical physics. These models have always proved to have a hidden 
algebraic structure which is responsible for their integrability [4]. The study of this 
algebraic structure also helps us to generalize the original model in many different ways. 
Motivated by the desire to understand the model (2) in this framework we have studied 
this system from a purely algebraic point of view. 

Our result is that when written in terms of suitable coordinates the Hamiltonian (2) 
takes a very simple form, which is nothing but the Hamiltonian of N Non-Interacting 
harmonic oscillators. All the results of Calogero [1] and Calogero and van Diejen [2] can 
then be derived without doing any calculations. Furthermore we will find the integrals of 
motion which are in involution, hence we prove the integrability of the system and find 
the explicit form of the action-angle variables . 



2 The New Dynamical Variables 



Define the variables 



aj = (^)*j vqi (10) 
a* = (p-)h e - iv v (11) 
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Since the parameter v has the dimension of inverse lenght, the above variables will have 
the dimension of (Action) 2 and will have the following poisson brackets: 



{aj, a k } = {a*, a* k } = {aj,a* k } = i5 jtk 



(12) 



The Hamiltonian will then take the following form: 

N-l N-l 

H = NXu a*cij + /iv a*dk 

j=Q j,k=0 

which can be put into the compact matrix form: 



(13) 



H = N\uA*A + fiuA^CA 



(14) 



where 



/ \ 

Oi 

a 2 



A = 



c = 



( \ 

1 1 ... 1 



1 1 ... 1 
1 1 ... 1 



1 1 ... 1 



(15) 



where dot stands for 1. Having all of its elements equal to 1, the matrix C has the 
following orthonormal eigenvectors: 



Cu a = £ a u a a = 0,1,2,. ...N-l 



(16) 



where 



Zo = N £i = & = £3 = .-6v-i = 
6 



(17) 



and 



(u a ) s = — != u as with uj = e 2 v 



N 



(18) 



i.e: 



u 



V 1 / 



/ \ 

1 



Ui = 



u) 



N 



( \ 
1 



U 2 



UJ 



UJ 



2N 



etc. 



(19) 



Clearly the matrix R which diagonalizes the hermitian matrix C is unitary, with 



RCR ] = D = diagonal{N, 0, 0, 0...0) 



(20) 



where the explicit form of R is given by 



Ra,f3 — (u Q )/3 — 



N 



UJ 



Defining the new variables B = RA , i.e: 



b„ = -4= cu as a s b* a = -^= u- as al 



N 



N 



(21) 



we find: 



H = N\vB ] B + jiuB^DB 



(22) 
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or 



H = N\u(b*bo + blh + b* 2 b 2 + .... + 6^_i&jv-i) + /i^^o (23) 

It is essential that the new variables satisfy canonical commutation relations which is a 
simple consequence of unitarity of R. 

{b a , b p } = {&• , b;} = o {b a , b;} = i8 a , p (24) 

Defining now the quantities I a = b* a b a which are global functions of the old coordiantes 
we find : 

{I a ,Ip} = Va,fl (25) 

The quantities I a are the N integrals of the motion, in involution with each other and 
the Hamiltonian is a function of these integrals, i.e: 

H = N\u(l + h + I 2 + ....Jjv_i) + tiuNI (26) 

Therefore the system (2) is integrable in the Lioville sense. The explicit form of the 
integrals I a is found from (10-11) and (21) to be : 

1 I 2vra 

4 = 7T7r Yl (PkPk') 2 cos[v(q k , - q k ) + - k)] (27) 

k,k'=0 

These functions , having the dimension of action are in fact the action variables. We 
can also find the angle variables Q a . These are: 

Q a = \xn-^ (28) 



It is straightforward to check from (24) that the action-angle variables have canonical 
poisson brackets: 

{la, If,} = {Qa, Q P ] = {Qa, Ip} = 8 a ,p (29) 

The total momentum P has a simple expression in terms of the action variables: 

P = is(l + I 1 + I 2 + ....I N - 1 ) (30) 



It is now clear that the initial value problem can be solved in very simple and closed 
form. The equations of motion are: 

^-b a = ivN\b a a^O (31) 
at 

^-bo = ivN(X + fi)bo (32) 
at 

with solutions 

b a (t) = b a (0)e mxt a^O (33) 

and 

b (t) = b (0)e iN{x+ ^ (34) 

Using the inverse transformation of (21) one obtains the time evolution of the varibles 
(aj,a*) and hence of the original coordiantes and momenta. It is only for the difference 
between evolution of b a ^o and b that the varibles (ai,a*) and hence qi and pi have a 
complicated-looking evolution. Otherwise ( in the case \i = ) one will find that the 
variables (a i: a*) have exactly the same simple time evolution as the variable (h,b*). 
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3 Generalizations 



One can now generalize the system (2) in many different ways. In fact from our analysis 
it is seen that if the matrix C in (14) is replaced by any other hermitian matrix , none of 
the main results of this paper will change, except the explicit form of the diagonalizing 
matrix. A good choice of the matrix C which is perhaps much better than (15) on 
physical grounds , is the following: 



C = 



1.. 
10 1. 
.10 1 



\ 



1 1 



V 



a . . 1 

where dot stands for 0. This matrix leads to the following Hamiltonian 



/ 



N-l N-l 

H = NXu a*cij + [w ^2 a *j a j+i + a* +1 cij+ 

j=0 j=0 



(35) 



(36) 



or 



N N-l 
3 3=0 

^(piPn)^ cos[u(q 1 - q N )} 



(37) 



which represents a system with local ( nearest neighbor ) interaction , rather the long 
range interactions implied by (2). Physically, systems with local interactions are much 
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more interesting than those with long-range interactions. Also it is meaningful to think 
of (35) as a one dimensional system whereas for the system (2) where all the particles 
interact with each other the concept of dimensionality is somehow vague. The parameter 
a in (34-36) is used to treat systems with various boundary conditions. When a = we 
are dealing with open boundary conditions and when it is non-zero we are dealing with 
periodic boundary conditions. 

4 Discussion 

We have shown that the integrable system proposed by Calogero when written in suitable 
coordinates is in fact a system of non-interacting harmonic oscillators. Our results also 
show how (in the quantum case) the spectrum of this model can be calculated in a quite 
simple way. From (24) it is seen that the energy eigenvalues are given by : 

E = NXu(n + rii + n 2 + ■■■■ n N-i) + Nnvn (38) 

where rij 's are positive integers. At first sight, it may seem disappointing that the 
solvable systems of [1,2] are in fact non-interacting harmonic oscilators in a different 
guise. However one should note that after all any integrable system when written in 
terms of the action angle variables will be a non-interacting system. Therefore from a 
rigorous point of view the systems studied in [1,2] deserve to be named " new solvable 
systems " although they are too easy to solve when treated properly. Other integrable 
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systems of the form 

N 

H = J2 (PjPk){^ + ^cos[u(qj - q k )}} 

j,k 

have been proposed by Calogero in [5]. The integrable structure of these systems have 
been analysed algebraically in [6,7] , where it has been shown that these systems are 
nothing but a system of spins with long range interaction. 
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